Threshold Phenomena under Photo Excitation of Spin-crossover Materials with 

Cooperativity due to Elastic Interactions 
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Photo-induced switching from the low-spin state to the high-spin state is studied in a model of 
spin-crossover materials, in which long-range interactions are induced by elastic distortions due to 
different molecular sizes the two spin states. At a threshold value of the light intensity we observe 
nonequilibrium critical behavior corresponding to a mean-field spinodal point. Finite-size scaling of 
the divergence of the relaxation time is revealed by analysis of kinetic Monte Carlo simulations. 

PACS numbers: 75.30.Wx 75.50.Xx 75.60.-d 64.60.-i 



Molecular solids whose molecules can exist in two 
electro-vibrational states, a low-spin ground state (LS) 
and a high-spin excited state (HS), are known as spin- 
crossover (SC) compounds. Due to higher degeneracy 
of the HS state, such materials can be brought into a 
majority HS state at low temperatures by several meth- 
ods, including light exposure, known as Light-Induced 
Excited Spin-State Trapping (LIESST) Q. If the inter- 
molecular interactions are sufficiently strong, this change 
of state can become a discontinuous phase transition such 
that the HS phase becomes metastable and hysteresis oc- 
curs. The different magnetic and optical properties of the 
two phases make such cooperative SC materials promis- 
ing candidates for switches, displays, and recording me- 
dia 0, and the cooperative properties of SC materials 
have therefore been studied extensively In addition 
to SC materials, models with similar long-range inter- 
actions are relevant to a range of physical phenomena, 
including earthquake faults [3j and phase transitions in 
polymers 

In a typical model of a SC material, the energy differ- 
ence between the HS and LS states is taken as D > 0, 
while the degeneracy <?hs of HS is larger than t/LS, such 
that HS is preferable at high temperature. We can ex- 
press the HS (LS) spin state of the i-th molecule by 
o~i = +1 (—1), giving the infra-molecule Hamiltonian, 
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where g = gus/ghs- Here the difference of the degen- 
eracies are expressed by a temperature dependent field. 
In order to provide the cooperative property in the SC 
transition, a short-range Ising-type interaction has pre- 
viously been adopted [a, Q • Recently, we pointed out an 
alternative mechanism with clastic interactions due to 



the volume difference between the LS and HS molecules 
0, [1] • The elastic interaction mediates the effects of lo- 
cal distortions over long distances, and it has been found 
that the critical behavior belongs to the mean-field uni- 
versality class |9|. 

Here, we study the dynamics under photo-excitation 
from the perfect LS state at low temperature. For a single 
molecule, we find a smooth increase of the HS fraction, 
and the saturated HS fraction is a smooth function of the 
irradiation intensity. On the other hand, if we include a 
cooperative interaction, a kind of threshold phenomenon 
appears. When the irradiation is weak, the system stays 
in the LS state. If the irradiation becomes stronger than 
a threshold value, the system jumps to the HS state. 
Thus, the stationary state suddenly changes from LS to 
HS at a threshold value of the irradiation intensity. This 
sudden change is analogous to the sudden change of mag- 
netization at the coercive field in the hysteresis loop of a 
ferromagnet following reversal of the applied field, known 
as a "spinodal point." 

The dynamics of field-driven first-order phase tran- 
sitions have previously been studied in detail for 
Ising/lattice-gas models with short-range interactions. In 
that case, the nucleation and growth of clusters play 
important roles, and different regimes have been iden- 
tified, depending on the relative time scales of nucle- 
ation and growth [Io| . If the average time between nu- 
cleation events is long compared with the time it takes 
a single cluster to grow to a size comparable with the 
system, the transformation will happen stochastically in 
a Poisson process via a single cluster (single-droplet, or 
SD, regime), and the average lifetime of the metastable 
phase is inversely proportional to the system volume. If 
the average time between nucleation events is short com- 
pared with the growth time, many clusters nucleate while 
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FIG. 1: (Color online) Time series of the order parameter 
M(t) under irradiation. Irradiation intensities a = 0.001, 
0.002, 0.003, 0.004, and 0.005 from bottom to top. We plot 
10 samples for each value of a. The system size is L = 50. 
The time t is measured in Monte Carlo steps per site (MCSS). 



the first one is growing, resulting in an almost deter- 
ministic process (multi-droplet, or MD, regime), and the 
metastable lifetime is independent of the system volume 
(Avrami's law) [Io| . 

In this paper, we study the photo-induced transition in 
a model of SC materials with elastic interactions, which 
we have developed previously d, 0]. In this model, we 
found that the appearance of large, compact clusters 
is strongly suppressed, and the system remains uniform 
near the critical point and also during simulations of field- 
driven hysteresis [s[ . In the process of the photo-induced 
transition, we also expect to see effects of the effective 
long-range elastic interactions. Here, we therefore inves- 
tigate the relaxation time near the threshold point and its 
dependence on the system parameters. We analyze the 
dependences from the viewpoint of the mean-field spin- 
odal, which shows qualitatively different properties from 
field-driven phase transformation in short-range models 
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Here we use Monte Carlo (MC) simulations according 
to the constant-pressure method In addition, we in- 
clude a process of photo-excitation by the following pro- 
cedure. We choose a site randomly and change it to the 
HS state with a probability a if it is previously in the LS 
state. We denote a as the irradiation intensity. 

The order parameter for the present model is the frac- 
tion of HS molecules, /hs> which is expressed by the 
"magnetization" 
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Here we adopt the two-dimensional model of Ref. with 
the same parameters, including g = 20 and D = 0.5. 

In Fig. [T] we depict typical time series of the HS frac- 
tion under irradiation. We find a sudden change in the 
saturated values. For a < 0.003 the system stays in the 
LS state, while it jumps to the HS state for a > 0.004. In 
Fig. [21 we show the dependence of the saturated value on 




FIG. 2: (Color online) Saturated order parameter M, propor- 
tional to the HF fraction, shown vs the irradiation strength 
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FIG. 3: (Color online) (a) An example of the time evolution 
of the order parameter M(t) under photo irradiation with a = 
0.004. Time is given in MCSS. (b) Spin configuration during 
the process of photo irradiation shown in (a) at t — 6000 
MCSS. 



a. Because of the cooperative interaction, the HS frac- 
tion remains low when the irradiation is sufficiently weak. 
On the other hand, it approaches unity when a exceeds 
a threshold value. This threshold value of a is estimated 
from Fig. [H to be slightly below 0.004. 

We also find a sudden change of M as a function of 
time. We depict an example of the relaxation process 
for L = 50 in Fig. [3]Ja), which shows a sharp change 
around t = 5000 Monte Carlo steps per site (MCSS). 
The HS fraction stays at an intermediate value (M ~ 
—0.2) for a while, and then it rapidly increases to the 
saturated value. This sudden change is reminiscent of 
a transformation via single-droplet nucleation from the 
metastable phase to the equilibrium phase. However, no 
visible clusters are observed. In Fig. [3^b) we depict a 
configuration at t = 6000 MCSS, near the starting point 
of the sudden change. This absence of cluster structures 
is the same characteristic which we found previously in 
simulations of hysteresis in this model Q . 

Next, we study the system-size dependence of the re- 
laxation time. Here we define the relaxation time r as the 
time at which the average of Cj reaches 0.5, correspond- 
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TABLE I: List of average and the standard deviation of the 
distribution of the relaxation time. 
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ing to a HS fraction of 0.75, starting from the perfect 
LS spin state (<Tj = —1). We perform -/V samp i e (> 200) 
independent runs for each parameter set and average the 
results. In Fig. [4] we depict the size dependences of the 
relaxation times r for a £ [0.00365, 0.00400], obtained by 
averaging over the samples. Here, we plot the data as a 
function of L 2 / 3 because it is expected that the relaxation 
time diverges as L 2 / 3 at the spinodal point (see below). 

We found the sudden change of the saturated value of 
M around asp ~ 0.0385 in Fig.O Below this value, the 
relaxation times increase rapidly with L, while above it 
the relaxation times saturate for large L. As shown in 
Fig. [TJ for a < asp the times of the jump are widely dis- 
tributed. This is a sign of stochastic nucleation, although 
we do not observe any clusters in Fig. [31(b) . For a > asp, 
the transformation process becomes nearly deterministic. 

It is also found that the size dependence of the stan- 
dard deviation of the relaxation time \J (r 2 ) — (r) 2 is 
very different below and above the threshold. Below asp 
the standard deviation is about the same as the mean 
value while it decreases with increasing L above asp- The 
former suggests a kind of Poisson process in the SD re- 
gion and the latter a kind of Gaussian process in the MD 
region. In Table I, we list typical examples. 

It should be noted that the relaxation time increases 
with the system size. This system-size dependence of 
the relaxation time is qualitatively different from the be- 
havior in systems with short-range interactions [13]. In 
the latter case, the relaxation time is determined by 
the probability of creating a critical localized, compact 
droplet, whose free energy is independent of the system 
size. Droplets thus nucleate in a Poisson process of rate 
L d I, where d is the spatial dimensionality of the system 
(here, d = 2). Thus, the relaxation time becomes short 
when the size increases in both the SD and MD regions. 

Compared to this dependence in short-range systems, 
the increase of the relaxation time in the present system 
with the system size is notable. The divergence when 
L — > oo indicates that the present threshold phenomenon 
is not just a crossover as in short-range systems, but a 
true critical phenomenon. Thus we will study its crit- 
ical properties by analyzing the size dependence of the 



FIG. 4: (Color online) Dependence of the relaxation time r 
on system size L and irradiation strength a. Shown as r 
vs L 2/3 for a = 0.00365, 0.00370, 0.00375, 0.00377, 0.00378, 
0.00379, 0.00380, 0.00381, 0.00382, 0.00383, 0.00384, 0.00385, 
0.00386, 0.00397, 0.00390, 0.00395 and 0.00400 from above to 
below. Straight-line behavior consistent with (JSJ) is seen for 
a = 0.00385 « a S p. 



relaxation behavior. 

The simple van Hove type dynamics using the mean- 
field approximation implies that the waiting time before 
a jump over the free-energy barrier for a < agp has the 
form 



exp(ciV), 



(3) 



where N = L d is the volume of the system and c is 
positive. The present case appears similar. Indeed, the 
present threshold behavior is essentially the same as the 
sudden phase change near a mean-field spinodal point, 
and we expect that the parameter dependence is the 
same as that near the spinodal point of the Ising ferro- 
magnet with weak long-range interactions (the Husimi- 
Temperley model), with the irradiation strength a here 
playing the role of the magnetic field. Thus, the relax- 
ation time for a > agp should diverge as 
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independent of L. At the spinodal point, the relaxation 
time diverges with the system size as 



(5) 



where N is the total number of spins [l4 , 15 , lit ]. 

We plot the relaxation times for various values of a 
as functions of N 1 ^ 3 — L 2 / 3 in Fig. 0] The relaxation 
times at a = 0.00385 are seen to follow a straight line, 
consistent with the properties of the mean-field spinodal 
point The relations ©, ||3J|, and §5§ can be combined 
into a scaling form, 
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FIG. 5: (Color online) Scaling plot of the scaling function 
for the relaxation time, In f(x) with x 3 ^ 2 . Here we adopt 
asp = 0.00385. The symbols denote the data of the same 
values of a as in Fig. [4] 
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FIG. 6: (Color online) Scaling plot of the ratio of the stan- 
dard deviation and the mean of the relaxation time with 
x = L 4,/3 (a — asp) with asp = 0.00385. The symbols denote 
the data of the same values of a as in Fig. [4] 



where the scaling function f(x) has the following asymp- 
totic dependences on the scaling variable x — L 4 ^ 3 (a — 
asp): 

{exp(|a;| 3 / 2 ) for x <C — 1 
x° for|x|«l • (7) 

x^ 1 ! 2 for i » 1 

In Fig. we plot the scaling function j{x) = t/L 2 ^ 3 on 
logarithmic scale vs |ir| 3 / 2 in order to have all the data 
including the large relaxation times at small values of a 
in one figure. We find that all the points collapse onto a 
scaling form, which shows the asymptotic behavior given 
in ([7]) . This scaling form has been found for the spinodal 
phenomenon in the Husimi-Temperley model as well [l6| . 

We also study the finite-size scaling of the standard de- 
viation of the distribution of the relaxation time. Above 
agp, the distribution of the relaxation time is narrowly 
localized, and the standard deviation decreases with in- 
creasing system size. On the other hand, far below asp, 
the distribution is exponential, so that the standard de- 
viation is of the same order as the average. In Fig. [5] we 
show the ratio of the standard deviation to the mean of 
the relaxation time vs the scaling variable x. The results 
collapse quite well onto a scaling function. 

In conclusion we note that the present spin-crossover 
model with elastic interactions displays a sharp spinodal 
point with critical properties consistent with those of the 
mean-field universality class. We hope that this truly 
critical nonequilibrium phenomenon soon will be experi- 
mentally observed during photo-irradiation of SC mate- 
rials. 
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